function dataPrep_impurityA_multi_filling()
%-1~0
label = {'-5/6','-19/24','-3/4','-17/24','-2/3','-15/24','-3/5','-7/12','-13/24','-1/2','-11/24','-5/12','-2/5','-9/24','-1/3','-7/24','-1/4','-5/24','-1/6','-1/12'};
nu = -[5/6, 19/24, 3/4, 17/24, 2/3, 15/24, 3/5, 7/12, 13/24, 1/2, 11/24, 5/12, 2/5, 9/24, 1/3, 7/24, 1/4, 5/24, 1/6, 1/12];
ind = [61,79,62,78,110,77,112,64,76,114,75,66,113,73,111,72,68,71,69,70];
lb = 25.66 / 6^0.5;
EC = 56.1 * 6^0.5;

%% calculate, plot, save data
data_to_save = [];

for i=1:length(nu)
    %record filling factor (value)
    data_to_save(i+1,1) = i-1;
    
    %calculate data
    result = CS1_0_extraction_2D_Phi(ind(i),0);
    
    data_to_save(1,4:length(result.az_r)+3) = result.az_r / lb;
    data_to_save(i+1,4:length(result.az_r)+3) = (result.az_phi - result.avg)*1000 / EC;
end

data_to_save = num2cell(data_to_save);
data_to_save(1,2) = {'label'};
data_to_save(1,3) = {'filling_factor'};
for i=1:length(nu)
    data_to_save(i+1,2) = label(i);
    data_to_save(i+1,3) = {nu(i)};
end

writecell(data_to_save,'processed data/multi_filling_2D_r.csv')


end
